Zero energy Majorana modes in spin ladders and a possible realization of the Kitaev 

model. 
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We show that in double-chain Mott insulators (ladders), disordered alternating ionic potentials 
may locally destroy coherence of magnetic excitations and lead to the appearance of spontaneously 
dimerized islands inside the Haldane spin-liquid phase. We argue that a boundary between the 
dimerized and Haldane phases of a spin-1/2 ladder supports a localized zero-energy Majorana 
fermion mode. Based on these findings we suggest a realization of a generalized Kitaev model 
where Majorana fermions can propagate in more than one dimension. 



INTRODUCTION 

A possible existence of Majorana zero modes in certain 
condensed matter systems has been a subject of much 
discussion in recent literature. The interest is caused by 
a strong nonlocality of Majorana zero modes which rises 
hopes of possible applications in such areas as quantum 
computation [1] and electron teleportation 2]. It has 
been argued that stable Majorana zero energy modes 
(MZEM) may appear on surfaces [3] and in vortex cores 
of certain superconductors with a nontrivial topological 
invariant (this includes spinless Px + ipy ones)[j|, Q. As 
far as magnetic systems are concerned, the existence of 
MZEM has been suggested only for a model system (the 
Kitaev model) [6| which material realization has not yet 
been achieved. 

This paper presents two main results. First, we show 
that localized MZEM appear on a boundary between 
spontaneously dimerized (SD) and magnetically disor- 
dered phases in systems of two coupled sgin-_l/2 chains 
(spin ladders). As was demonstrated in [7], 8], the SD 
phase may exist in spin ladders with sufficiently str ong 
four-spin interaction. There are reasons to believe [10| 
that such interactions are present, for instance, in spin 
ladder material CaCu203. Second, we suggest that a bal- 
ance between competing spin-disordered and SD phases 
can be shifted by application of a staggered ionic po- 
tential which effectively enhances the four-spin interac- 
tion. This may be important for applications since real- 
ization of such potentials is quite feasible by chemistry 
means. Based on these findings we suggest a particular 
realization of a generalized Kitaev model where MZEM 
are located on a manifold of arbitrary dimension. As in 
the original Kitaev model, these dispersionless modes fa- 
cilitate propagation of Majorana fermions with nonzero 
dispersion so that such propagation indeed becomes be- 
comes akin to teleportation as described in 

The idea that an alternating ionic potential can lead to 
SD was expressed in [ll|,[l^ where the authors studied 



insulator (IBI) in one dimension. It was found that since 
the Mott and IBI states differ under a local parity inver- 
sion, the transition driven by a change in the strength of 
the alternating ionic potential proceeds through an inter- 
mediate spin liquid phase characterized by a spontaneous 
dimerization. In the present paper we generalize these 
ideas for the case of (i) disordered ionic potentials, (ii) 
the systems consisting of two half filled chains coupled 
together. In the second part of the paper we construct a 
generalized Kitaev model. 



STAGGERED IONIC POTENTIAL IN A SINGLE 
HUBBARD CHAIN 

Let us briefly recall the single Hubbard chain case. At 
half filling such chain has a spectral gap in the charge 
sector and gapless excitations in the spin sector. At low 
energies the spin sector decouples from the charge one 
and can be considered separetly. For reasons of conve- 
nience we will use the continuous approximation for the 
entire Hubbard model. This is legitimate if the Mott- 
Hubbard gap is much smaller than the bandwidth. The 
effective Euclidian Lagrangian of the Hubbard model in 
this case is 
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a transition from the Mott to a conventional ionic band 



where $c.s are charge and spin bosonic fields, Vc.s are 
their velocities, K is the Luttinger parameter and G is 
the Umklapp coupling constant. All these parameters 
are related to the bare parameters of the Hubbard model. 
The cosine term in ([2]) is relevant for i^T < 1 (repulsion) 
and opens a gap (Mott- Hubbard gap Amh) in the charge 
sector. The spin sector remains gapless and the correla- 
tion functions of the spin fields display power-law decay 
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at T = 0. We omitted in ([3]) a marginally irrelevant term 
contributed by backscattering of the electrons. 

Let us consider the influence of a staggered ionic po- 
tential. Using the well known textbook bosonization for- 
mulae 13|, 14| for the staggered charge density we ob- 
tain the following contribution of such potential to the 
Lagrangian density: 



J dx F(^,a;)p3tag(x) (4) 
^ J dx y(7r, a;) sin -\/27r$c cos v^^s, 

where V{tt, x) is an envelope made of the Fourier harmon- 
ics of the potential concentrated around wave vector tt, 
and yOstag (x) is the continuum limit part of the staggered 
density (— l)"C+rC„o-- Since the charge sector has a gap 
and the vacuum average of sm\/2n^c is zero, the effect 
of the staggered potential shows up only in the second 
order of perturbation theory in V{2kF = tt) |11|]: 



(5) 



6H = X d^^sRA^sL - {2iTa^y cos V87r$ 



where ^sR;L are chiral components of the scalar field <i?s, 
\ ^ V'^ and a is a short-distance cutoff of the bosonic 
theory. More precisely, the coupling A is expressed as a 
convolution of the disorder average of the potentials with 
the correlation function of the charge fields: 



A(a:) j drdy F(7r, X -I- y)V^(7r, y) X (6) 
((sin %/2^$c(t, x + y) sin %/2^$,(0, y))) (r^ + y'^/vl), 

and in principle is coordinate dependent. To provide the 
convergence of ([6]) in time domain we need an exponential 
decay of the dynamical correlation function which occurs 
only if the charge gap is nonzero. To estimate the average 
coupling we may assume that 



(F(7r,x)T/(^,0)) =F2exp(-|a;|/0. 



(7) 



In the most interesting case when the disorder correlation 

length ^ ^ Vc/ ^ we have A ~ ^a^/\l ■ Interaction ^ is 
marginally relevant and competes with the marginally 
irrelevant interaction of the same form present in the 
spin sector of the original Hubbard model. Therefore 
such interaction will open a spin gap only when A ex- 
ceeds some critical value related to the aforementioned 
marginally irrelevant backscattering term. Refs. 12| 



consider only a uniform potential, and one could get an 
impression that the onset of the gapped dimerized phase 
would require doubling of the unit cell across the whole 
chain. However, from the above example it follows that, 
to achieve the desired effect, it is sufficient to break trans- 
lational symmetry only locally. 



DOUBLE CHAINS. THE CASE WHEN CHARGE 
GAPS >> THAN THE SPIN ONES 

Now we consider a system of two coupled half-filled 
chains. For simplicity one may think about them as Hub- 
bard chains even though precise details of the charge sec- 
tor are not important. It was shown in that at half 
filling, when the low-energy dynamics in the spin sector 
of a single chain coincides with the one for the spin S=l/2 
Heisenberg model, the continuum limit description of the 
ladder is given by the model of four Majorana fcrmions 
with the Lagrangian density 



a=l 



(8) 



where the Majorana fields xRj Xl in the path integral 
are real Grassmann numbers and parameters v,mt,ms 
depend on the bare interactions. This description works 
well when the interchain exchange is much smaller than 
the exchange along the ladder. Then parameters mt,ms 
(masses of the triplet and the singlet Majorana fermions) 
are proportional to the linear combinations of the con- 
ventional J± and the four-spin Jcycie exchange integrals 
i,i: 



cycle; 



cycle- 



(9) 



{A,B are numerical coefficients). The spectra of the 
triplet and singlet Majorana fermions are et,s{k) = 

^ (ufc)2 -|- mf they do not depend on the signs of mt^s- 

However, since the spin operators are nonlocal in terms 
of the fermions, their correlation functions do depend on 
these signs. In particular, if rritms < the ladder is in 
the spin-disordered (Haldane) state characterized by co- 
herent magnetic excitations with spectral gap \mt\. If, 
on the other hand, rritms > the ground state is spon- 
taneously dimerized (SD) , that is the energy density ac- 
quires a nonzero staggered component. Away from quan- 
tum critical points where some of the masses are zero, the 
spin excitations have spectral gaps but become incoher- 
ent Q . In most experimental realizations of spin ladders 
the interchain exchange is either equal or larger than the 
exchange along the ladder. In that case the singlet gap is 
large and the corresponding excitation has no influence 
on the physical properties. However, in such systems as 
CaCu203 where the rung exchange is much smaller than 
the leg one the situation is different and the singlet mode 
(though appearing only in the magnetic continua) visibly 
affects the dynamical correlation functions (lo| . 

One may expect that by arranging a proper coordinate 
dependence of the interchain interactions one can create a 
situation when some or all of the Majorana masses change 
sign. This will lead to creation of Majorana zero modes 
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|15| . In [16j it was shown that this does happen when 
one puts a single vacancy on one of the chains; then all 
Majorana masses change sign. The result of such change, 
as it might be expected, is a creation of local spin 1/2 
localized at the vacancy. Here we would like to consider 
a different possibility, namely, when only an odd number 
of Majorana modes change sign. This happens either on 
a border between SD and magnetically disordered state 
or, if the SU(2) symmetry in the spin sector is lost, by 
a change of sign in one the components of the Majorana 
triplet. The latter can be achieved by applying an easy 
axis anisotropy. 

To be definite we will consider the SD case. We suggest 
that from the practical point of view a local change of to^ 
can be achieved by applying staggered ionic potentials 
acting on the staggered charge densities of both chains: 



(10) 



1=1.2 



where Vi{x) is a slow function of a; = naQ. In order to 
achieve the desired effect that only rris changes sign we 
need V1V2 < 0- In the opposite case V1V2 > the sign 
change will occur in the triplet mass mt. For simplic- 
ity we consider only the former possibiity. To perform 
actual calculations we adopt the following hierarchy of 
energy scales 4i|| :s> V, Amh ^ t±, where Amh is the 
Mott-Hubbard gap for the individual chains, V is the in- 
terchain density-density interaction and t±_ is the inter- 
chain tunneling. With such assumptions the high energy 
degrees of freedom are the charge modes which can be 
described by the following model: 



(11) 



Vdx^cidx^c2 - G {cob VStt'I'cI + cos \/87r$c2 



where fields $ci,'I*c2 are the charge fields of individual 
chains, V being proportional to the interchain density- 
density forward scattering. Passing in ([TU)) to the con- 
tinuum limit and bosonizing the resulting expression we 
obtain: 



5H - Fi(a;)sin(V2^$ci)(e4 
-I- V2(x)sin(%/2^$c2)(e4 



£-) 



(12) 



where e± are the symmetric and antisymmetric combi- 
nations of the chain spin-dimcrization operators [l7| |. In 
the Majorana representation (jS]) these operators are ex- 
pressed as linear combinations of products of four order 
or disorder parameters of the corresponding Ising models 
(each noninteracting Majorana fermion is equivalent to 
the quantum Ising model; see 01 and also 13 1 for details): 



At small distances the following fusion rules take place: 
e±(R-^p/2)e±(R-p/2)= (14) 

i(?)n[^^^^'° ^.(R)+---] 

i=i 

= const ± —nra >^ ir^ap > + ■ ■ • , 

^ i=i ^ \i=i / 

where — X^rXl- Using fusion rule (1141) and keeping 
only the most relevant terms we obtain in the second 
order in the potentials the following shift of all masses: 



ms^t = mfl + Ci2{x), 



Ci2{x)^ / AydrViix + y)V2{v) x 



(15) 



((sin V27r$ci(r, x + y) sin \/27r$c2(0, y))). 

Naturally, C12 vanishes in the absence of interchain in- 
teractions. Since in the spin-liquid state the masses 
of triplet and singlet Majorana fermions have different 
signs, such shift may lead to a change in sign of one of 
the masses. In the uniform case this would push the sys- 
tem to a spontaneously dimerized phase via Quantum 
Critical Point [8]]. If the ionic potential is coordinate 
dependent the sign change will occur locally. For the 
Hubbard model where the interchain coupling is antifer- 
romagnetic only the singlet mass changes sign. Then on 
phase boundaries between SD and disordered phases a 
single Majorana mode has zero energy states, as it occurs 
with fermionic states located on the boundary of quasi- 
ID p-wave superconductor Q . Another direct analogy is 
a boundary between ordered and disordered states in the 
quantum Ising model. 

A purely one-dimensional case (single ladder) will dis- 
play disorder effects studied in 18 1 and [31 where the 
case of random potential was considered. It was shown 
that when such boundaries are randomly distributed with 
a final concentration rii they give a singular contrubution 
to the specific heat: 



C{T)^ndn-\TolT). 



(16) 



MlM2A'3M0, 



(TlCr20'3Cro- 



(13) 



We remark that the models considered in [18|,[16| dealt 
with the situation when the number of Majorana zero 
modes on each domain wall was four. In that case, as we 
have already mentioned, each domain wall carries local 
spin 1/2. 



HOW TO CONSTRUCT THE KITAEV MODEL 

A much more interesting situation emerges when the 
number of modes is odd like in the case we consider. Let 
us consider the above ladder model with large areas of the 
dimerized phase separated by large areas of the Haldane 
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FIG. 1: Construction of the Kitaev model. Double rails rep- 
resent spin ladders, the ellipses denote interactions (I18|l facil- 
itated by the magnetically active media of the substrate de- 
scribed in the main text. Areas with different filling represent 
SD and Haldane phases. MZEMs are located at boundaries 
between the phases. 



ixlxl + XRXR)y.xAxlx\ + XRXR)y+l,x^, 

where coordinates mark positions of domain walls. For 
each wall one has to specify whether it holds a zero mode 
with right or left chirality and replace x° operator with 
its zero mode component: 



x°r{x) \ _ i{xc,y) ± ms(x')dx'/-« 

xl{x)) - U 



1 



.(20) 



phase, such that the wave functions of the zero modes 
from different domain walls do not overlap. Now let us 
put these ladders on top of a superparamagnet whose 
magnetic susceptibility has a peak at (tt, 0) wave vector. 
So the strongest magnetic fluctuations in such a material 
are smooth along the ladder direction and oscillate in 
the direction perpendicular to the ladders. Then they 
strongly couple to the 



K'-^-^ixlx'k + xlxll (a = 1,2,3), 



(17) 



the operator representing the difference between magne- 
tization densities of two legs of each ladder Q . Integrat- 
ing over the bulk magnetic fluctuations we obtain the 
effective interaction between Majorana modes from dif- 
ferent ladders: 



V^r 



x)Jab{x - x', r)K\y + r, x'), (18) 



Jab{x - x'- r) = {{N'^iy, x)N\y + r, x')}}, 

where y is a location of the ladder in the transverse direc- 
tion, X, x' are the coordinates along the ladder and N°- 
is the a-th component of the staggered magnetization of 
the paramagnet. Since the paramagnet is disordered, this 
interaction is short ranged. To get the Kitaev model we 
have to arrange the positions of domain walls in such a 
pattern that a wall on one ladder interacts only with one 
other wall. Such pattern does not need to be regular 
though disorder will lead to localization of the propagat- 
ing Majorana modes. 

As an example we may consider a pattern where each 
ladder contains alternating SD and Haldane phases of al- 
ternating lengths a and b. The ladders are shifted by 
amount c < a in the y direction (see Fig.l). All lengths 
are assumed to be much larger than the domain wall 
width so that the wave functions of MZEMs on different 
domain walls do not overlap. Since the correlation func- 
tion Jab is short range, we can neglect all interactions 
between domain walls except between the ones located 
on neighboring ladders. Hence we get 



H = 

E 

j/,a=l,2,3 



dx 



-i-x'kd^x'k + xld.xD+iratx'kxl 



where the sign depends on the sign of Smg — ms(oo) — 
r7is(— cxd), "f(xc,y) is a zero energy Majorana fermion 
mode satisfying the Clifford algebra anticommutation re- 
lations, Af is the normalization factor, the dots stand for 
contributions of higher energy modes. 

Model is a particular type of the Kitaev model 
fermionized by means of the Jordan- Wigner transforma- 
tion, as in [l§|. The latter method, when the original spin 
model is treated as a model of coupled chains and Jordan- 
Wigner transformation is done for each chain, represents 
certain advantages since it allows to avoid a cumbersome 
gauge fixing used in the original paper [1] . The most im- 
portant feature of the arrangement depicted on Fig. 1 
is that each MZEM in Hamiltonian (|19p interacts only 
with one nearest neighbor exactly as in the original Ki- 
taev model. In this way the products of neighboring zero 
Majorana modes ^{x,y)^{x,y + 1) are integrals of mo- 
tion and can be replaced by constants. Thus the four 
fermion term in (|19p becomes effectively a hopping term 
for Majorana fermions which now can propagate in 
the transverse direction. The latter fermions have spec- 
tral gaps, but such sector also exists in the Kitaev model. 

To make sure that MZEMs are protected, we consider 
the operators which couple these modes to local pertur- 
bations. The dangerous ones are those which can couple 
MZEMs located at different positions. Neither K op- 
erator considered above nor the staggered components 
of magnetization can do this since they contain massive 
modes and their correlations falls off exponentially. This 
leaves the staggered energy density operators e± which 
are products of four order ( disorder) Ising model oper- 
ators (|13p . The antisymmetric staggered energy density 
e„ contains operators Ua (a=l,2,3) which have no ground 
state average and therefore vanishes. Operator e+ has a 
nonzero ground state average affected by MZEMs and 
therefore may couple to it. The coupling originates from 
the change in /zq; the operators from the triplet sector 
do not experience any change at the boundary. Using 
Eqs. (83,84) from [11] we get (a; < is in the SD phase): 



(19) 



s+(x)) = (a;lm?m,)i/«exp{-/K(x-Xe)]},(21) 



f{x) = -e{x)x + - Ko{\x\) + K.^{\x\) 
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FIG. 2: e_|_(a;)/e( — cxd) as a function of 2ms (x — 



This function is depicted on Fig. 2. Here to^ stands 
for the asymptotic value of the singlet gap far from the 
boundary. As we see, at the boundary where the MZEM 
is located the symmetric staggered energy density van- 
ishes: (e+(a;)) ^ \x — Xc\^/^. This fact renders the inter- 
action with MZEM numerically weak. 



CONCLUSIONS 

We have demonstrated that zero energy Majorana 
modes (MZEMs) robust against external perturbations 
exist on a boundary between the Haldane and SD phases 
of spin S=l/2 ladders. We have also suggested a practical 
way to manufacture systems with such phase boundaries. 
Namely, the effective interactions in spin ladders can 
be manipulated by purely chemical means through stag- 
gered ionic potentials. There are two ways MZEMs from 
different domain walls interact between each other: by 
direct overlap and via interaction of a particular Fourier 
component of the magnetization (|17l) . The former inter- 
action leads to formation of an impurity band inside of 
the spin gap with associated singularity in the specific 
heat (|16p. However, by maintaining a proper distance 
between the phase boundaries one can make this overlap 
small. For this case we have suggested a particular ar- 
rangement when the interactions between MZEMs from 
different ladders create a network similar to the one ex- 
isting in the Kitaev model in its gapped phase. 
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